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Summary - 

This paper deals with the initiation, growth, and possible 
arrest of fracture in cylindrical shells containing initial 
defects. For those defects which may be approximated by a part- 
through semi-elliptic surface crack which is sufficiently shallow 
so that part of the net ligament in the plane of the crack is 
still elastic, the existing flat plate solution is modified to 
take into account the shell curvature effect as well as the 
effect of the thickness and the small scale plastic deformations. 
The problem of large defects is then considered under the assump- 
tions that the defect may be approximated by a relatively deep 
meridional part-through surface crack and the net ligament 
through the shell wall is fully yielded. . The results given in 
the paper are based on an 8th order bending theory of shallow 
shells using a conventional plastic strip- model to account for 
the plastic deformations around the crack border. 

For shells with shallow surface cracks the plasticity cor- 
rected stress intensity factor given in Section 2 of this paper 
may be used to study the phenomena such as plane strain fracture, 
fatigue crack propagation, and stress corrosion cracking. To 
study the fracture problem in shells containing large defects m 
Section 3 the results of COD calculations _ are presented. These 
results may be used in conjunction with either a critical crack 
openinq stretch type fracture criterion, or a plastic instability 
criterion which is described in Section 4. This section also 
’ contains a comparison of the calculated results with those 
obtained from the experiments on steel pipes. 

In Section 5 a preliminary attempt is made to come to grip 
with the important problem of dynamics of crack propagation in 
pressurized shells which is coupled with the problem of fluid . 
mechanics. As an example, the initial stage of crack propagation 
in a finite volume container is considered. The results indicate 
that in some cases the fracture arrest may be possible due to 
depressurization . 
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List of Symbols 


A : Area of crack opening (orifice) . 

A r A.: Membrane and bending components of the stress 
m intensity factor ratio in shells. 

2a : Crack length. 

COD: ... Crack opening displacement. 

d :* Depth of the surface crack. 

= 4aa Y /E : normalizing factor for 6. 

E : Young ' s modulus . 

h : Shell thickness. 

K : Stress intensity factors. 

N : Membrane load in circumferential direction, 

o 

p ' : Pressure. 

R : Shell radius of curvature. 

6 : Crack opening displacement or stretch. 

X = ' [12 (1-v 2 ) ] ly/4 a/ (Rh) •: shell parameter. 

Oy : Yield strength. 

o’ : Flow stress. - 


1. Introduction 

Generally, in pipes, pressure vessels, and a great variety 
of other plate and shell structures the main cause of structural 
failure appears to be the nucleation and the growth of cracks 
from an existing dominant imperfection or flaw which is usually 
a surface or an internal defect of arbitrary shape and orienta- 
tion. In some cases inspection and nondestructive testing pro- 
cedures may provide a fairly good estimate of the size, shape, 
orientation, and location of the dominant defect. However, in 
most analysis relating to the reliability and integrity of the 
structure, generally one simply conjectures the existence of such 
defects and proceeds with the application of a relevant fracture 
criterion. Similar to all failure studies, the work leading to 
the application of the fracture criterion has three . different 
components which may initially be considered independently. 

These are: 

(a) . The experimental determination of the "characteristic 
strength parameter" representing the inherent resistance of the 
material to fracture corresponding to the conjectured mode of 
fracture. 

(b) The solution of the mechanics problem for the given 
geometry and the external loads in order to determine the "criti- 
cal load factor" representing the severity of the applied loads. 

(c) Development of an appropriate fracture theory and a 
related quantitative "fracture criterion". This is usually 
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nothing but the direct comparison of the "characteristic strength 
parameter" of the material with the critical load factor^. Xt 
is clear that in fracture studies the development of the fracture 
theory and selection of a proper fracture criterion must precede 
the experimental determination of the characteristic strength 
parameters and the theoretical evaluation of the corresponding 
critical load factors. 

In shell structures, including the pipes and pressure ves- 
sels, if the flaw geometry and the external loads are such that 
the stress State in the medium is basically elastic and the 
plastic and other inelastic deformations are restricted to the 
relatively small toroidal region enclosing the crack periphery, 
the ensuing fracture is generally of plane strain type and may 
be classified as "low energy" type. Consequently, the "fracture 
toughness" G IC (or the corresponding "critical stress intensity 
factor" K xr ) , the fatigue crack propagation constants, the cor- 
rosion cracking threshold value of the stress intensity factor 
K TQrp , and similar quantities may be considered as the material 
parameters representing the resistance of the structure to fail 
ure. In this case the appropriate "load factor" is the stress 
intensity factor K (or quantities which may readily be evaluated 
in terms of K) , and hence in fracture studies it is sufficient 
to perform a linear analysis.' 

On the other hand, if the crack geometry, the external loads, 
and the environmental conditions are such that there are large 
scale plastic deformations around the crack front, the resulting 
fracture will be of "high energy" type. In this case one rele- 
vant criterion which lends itself reasonably well to standardiza- 
tion and analytical treatment is that of the "critical crack 
opening stretch" S r . The criterion appears to work rather suc- 
cessfully for "plane stress" type fracture of thin sheets with 
relatively large plastic zones around the crack. 

The calculated results for the crack opening displacement 
(COD) in shells indicate that for a given shell geometry COD is 
a monotonically increasing function of the hoop stress oqq with 
a slope which also increases quite rapidly for increasing values 
of aAfi/dvr Ov being the yield strength. This suggests that 
around a certain value of the hoop stress, any small increase in 
the applied load may cause a considerable increase in COD, result 
ing in, what may be called, "the plastic instability" of the 
material ahead of the crack and causing rupture. .Thus, the basis 
of another fracture criterion to be used in studying -the rupture 
of thin sheet structures may be the concept of plastic instabil- 
ity which may quantitatively be characterized in terms of a 
standard slope in the COD vs Cqq relationship. 


The primary objective of this paper is to provide the neces- 
sary analytical tools for the application of different types of 
fracture criteria mentioned above to the shell structures. Of 
particular interest will be the problem of fracture propagation 
and arrest in a cylindrical shell containing an initial surface 
flaw. It will be assumed that the flaw has the most unfavorable 
geometry and orientation, namely it is a crack lying on a meridi- 
onal plane. First we will consider the problem of a part- through 
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crack for which either the crack depth-to-shell thickness ratio 
d/h is sufficiently small or the external loads (and/or the tem- 
perature) are sufficiently low so that the plastic deformations 
are confined to the neighborhood of the crack border and the 
stress state in the shell is basically elastic. The main analyt- 
ical problem here is the evaluation of the (plasticity-corrected) 
stress intensity actor along the crack front. 

We will then consider the problem of a relatively deep sur- 
face crack or a relatively large internal crack for which the 
net ligament in the plane of the crack is fully yielded (Figures 
lb and 2) . In this case it will be assumed that the plastic 
deformations are also spread around the crack ends through the 
entire wall thickness. Here the main analytical problem is the 
evaluation of the crack opening stretch at the locations of 
expected fracture propagation. Calculated COD values will 
compared with the experimental results obtained by Fearnehough 
and Watkins 1 , and by Cowan and Langford . The results of the 
burst tests performed by Anderson and Sullivan , will be used to 
test the validity of the critical crack opening stretch as the 
fracture criterion. To study the effectiveness of the concept 
of "pLstic instability" as a correlation tool in "plane stress- 
type fracture, the calculated results will be compared with the 
experimental results obtained by Kiefner, et al. Finally, the 
problem <6f the rupture of the net ligament and the ensuing crack 
propagation and arrest in a pressurized cylinder of finite volume 
will be considered theoretically. 


2. Elastic Shell with a Part-Through Crack 

The elasticity solution of the part-through crack problem 
in plates or shells is not yet available. For plates there are 
some approximate solutions which are believed to give reasonably 
accurate results under some simplifying assumptions regarding 
the crack geometry and the boundary layer effect. These solu- 
tions (which are reviewed in the ASME volume on The Surface 
Crack s ) are essentially based on the modification of the elliptic 
crack problem considered by Green and Sneddon 6 . A brief review 
of these results with further modification taking into account 
the shell curvature effect may be found in the paper by Erdogan 
and Ratwani 7 . In this paper only the expression for the stress 
intensity factor in a cylindrical shell will be repeated. 

Consider the crack geometry shown in Figure lb. Let us ^ 
assume that the crack depth- to- thickness ratio d/h, the applied 
load level, and the ambient temperature are. such that the mam 
bulk of the shell cross-section is elastically deformed, and the 
plastic deformations are restricted to a- relatively small region 
along the crack border. The (plasticity-corrected) stress inten- 
sity factor along the crack front may then be approximated by 


K ($) 



c.c 
t c 

Q 



2 , , ■ 2 , » 1/4 

cos cp + sin <p) 


( 1 ) 


where 
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<b is the angular parameter measured from the x axis (see 
Figure 2) which defines the leading edge of the crack 
through the following parametric equations: 

x *= a cos 4 > , z = -d sin<£ , 

N /h is the hoop stress acting on the shell outside the 
® perturbation zone of the crack, 

d,a,h are the dimensions shown in Figure la, 

' u 1 * • ' ‘ 

.2 


N' 

,2 - - o 


Q » [E^ (k) - 0.212 


.22 
h o y 


1/2 

3 


(3) 


E(k) , (k 2 = 1 “ d 2 /a 2 ) is the complete elliptic integral of 

the second kind, 

cjy is the yield strength of the material, 

2 


1 + 0.122(1 - d/2a) ^ 

c - — 4 I/? ' K } * 

t Y-.v [1 - (d/a) 4 ] 17,4 
c c = <A m“ 1) H + 1 ' - 

— 1 . 2 5 X 

A = 0.481A + 0.614 + 0 . 386e * , 

m 

A = [12(1-v 2 )] 1/4 V 72 ' 

(Rh) W 


(4) 

(5) 

( 6 ) 
(7) 


R being the (mean) radius of curvature of the cylindrical 
shell. 

In this analysis Q accounts for the plasticity correction, c fc is 
the thickness correction factor, c c is the shell curvature 
correction factor with A m being the stress intensity factor ratio 
in the shell with a through crack. 

3 . p ar t-Through Crack with Fully-Yielded Net Ligament 

In this section a brief description of _ the results _ regarding 
the part- through surface crack in a cylindrical shell with fully 
yielded net ligament will be presented. As shown in Figure 2a 
the geometry of the internal defect is generally quite irregular, 
and in a given "defect zone” there may be more than one internal 
crack. Here we will assume, that the total defect zone may be 
replaced by a single part-through crack of equivalent size and 
orientation, and of idealized shape (Figure 2b). Further, it 
will be assumed that the net ligaments connecting the idealized 
crack to the inner and the outer surfaces of the shell are fully 
yielded and carry only membrane stresses of magnitude ffy/ the 
yield strength of the material (Figure 2c). Finally, it will be 
assumed that the plastic deformations have be ^. 

ends of the crack up to a distance p (Figure 2b) and the plastic 
region may be approximated by a plastic strip similar to tha m 
the Barenblatt-Dugdale model. The plastic strip xn the yield 
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zone a< x < a+p = a p is assumed to carry a membrane load Nyy - N 
and a moment resultant Myy = M where N and M are as yet unknown 
and x,y,2 are*the rectangular coordinates in axial/ tangential, 
and radial directions. The main objective of the analysis is 
the determination of the crack opening stretch along the crack 
boundary. 

The method of solution and extensive computed results 
regarding the stress intensity factor and the COD in cylindrical 
shells were presented in the previous work of Erdogan and 
Ratwani 8 9 9 ' 7 . In this paper we will repeat some of the relevant 
results from the previous work and present some new results on 
the application of the plastic instability concept. 

Calculated results for the crack opening stretch 6 are 
shown in Figures 3-9. Figure 3 shows the crack opening displace- 
ment 6 0 in the neutral surface of the shell at the mid-section 
of the crack (x=0, z=0) for a through crack. In this and in the 
subsequent figures the normalizing factor for COD is defined by 

d x f 4aa Y /E , (8) 

where E is the Young 1 s modulus. Figure 4 shows the conventional 
. crack opening stretch 6 evaluated in the neutral surface of the 
shell at the, (actual) crack tip (x=a, z=0) . Figures 3 and 4 
clearly indicate the influence of the shell curvature (as repre- 
sented by the shell parameter X defined by equation (7)) on the 
crack opening displacement. In Figures 5-8 6 is the crack open- 
ing stretch (which includes curvature and bulging as well as the 
membrane effects) at the leading edge of the crack in the mid- 
section (x = 0, z = h/2-d) for a cylindrical shell with an exter- 
nal part-through crack of depth d. The figures essentially show 
the load carrying capacity of the shell for a constant crack 
opening stretch 6^ = 6c • In these figures the small circles 
shown for d=h for various values of 6/di correspond to the load 
levels No/hcy giving these COD values at the crack tip x- a in 
the shell having a through crack of the same length 2a. Gener- 
ally the extrapolated values of the load N 0 /hay fo r constant 6 
curves at d=h are lower than the load corresponding to the same 
COD at the crack tip in the through crack case. This means that, 
the critical COD criterion is valid, at least theoretically 
,i*t is possible to have leak before burst (or before any addi 
tional crack extension) . 

The load carrying capacity for the shell with a through 
crack is shown in Figure 9 • Here X represents the crack length 
( see equation (7) ) , and 6 is calculated at x — a. 

Figures' io and 11 show a comparison of experimentally meas- 
ured COD values in zircaloy pressure tubes used in steam gener- 
ating heavy water reactors (Fearnehough and Watkins , and Cowan 
and Langford 2 ) with those calculated by using the same dimensions 
and the same material constants (see Erdogan and Ratwani for 
details) . Generally the agreement appears to be quite satisfac- 
tory. Figure 12 shows the correlation of the results of the 
burst tests performed by Anderson and Sullivan on 2014-T6 alumi- 
num cylinders. In these figures the open circles correspond to 
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( 9 ) 


the elastic stress intensity factor in the shell given by 

1/2 

K = (A m + V (p o R/h) (ira c 5 

where Am and Ab, respectively, are the membrane and bending com- 
ponents of the stress intensity factor ratio in a cylindrical 
shell with a through crack, p D is the burst pressure and ac is 
the value of & at burst. K p *, is the plasticity corrected K 
obtained from (9) by replacing a c by a c + p, p being the plastic 
zone s i ze , computed from the plastic strip model. The figure 
indicates that, for the material under consideration, compared 
to K= constant, or Kp& = constant, 6 = <S C = constant appears to be 
a better rupture criterion. 


4 . Fracture Criterion Based on Plastic Instability 


From Figures 3 and 4 and from the similar results given for 
the part- through crack in a cylindrical shell by Erdogan and 
Ratwani 7 it is easy to observe that for a shell with a given 
crack geometry in the neighborhood of a certain value of the 
applied load N 0 /hay any small increase in N 0 may cause a very 
large increase in the COD (6 or fit) . This simply means that 
about this particular applied load. the phenomenon which takes . 
place at the leading edge of the crack is similar to the necking 
phenomenon observed in a ductile tensile bar, that is, the 
material is subjected to plastic instability. Moreover, the 
slopes of 6/di vs N 0 /hcy curves increase so rapidly with increas- 
ing N ft that the load corresponding to the plastic instability 
may be determined from these figures within an acceptable degree 
of accuracy by assuming a hypothetically selected high slope 
(such as 20) • Once a critical (standard) slope (in 5 vs N o 
curves) is selected, the load carrying capacity of the shell (N 0 
vs X) may then be obtained from Figure 4 for the through crack 
and from the similar figures for the part- through crack. These 
theoretical results are shown in Figure 13 where a slope of 20 
in 6/di vs No/hcry curves was used to determine the plastic insta- 
bility load N 0 for a given X and d/h^_ In this and income of 
the subsequent figures the notation a rather than a y is used. . 

As suggested by Kiefner, et al.\ in this type of analysis it is 
more appropriate to use a "flow stress". o rather than the stand- 
ard 0.2% offset yield strength_0y, 0 being somewhat higher than 
ov. 0 may be selected as (a) 0 = (l+cOoy, 0 being a fixed per- 
centage, or (b) a = (0y + 0 u )/2, a u being the ultimate strength, 
or (c) a ■ oy+Auy, Aoy being a constant (taken as 10,000 psi by 
Kiefner, et al." for pipe steels they tested). , Needless to say, 
the theoretical results given by Figure 13 are not affected by 
any particular selection of a. 

Figures 14 and 15 show the comparison of the theoretical 
and experimental results. The experiments were performed by 
Kiefner, et al." on steel pipes (X52, X60V, X65C, X60C contain- 
ing a through crack. In these figures the experimental results 
were normalized by selecting 0 = a Y +10,000 psi _ as suggested oy 
the Authors . A comparison of some of the experimental results 
given by Kiefner, et al. 4 for the part-through crack in steel 
pipes (X52, X60C, X60V) and the theoretical results shown m 
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Fiqure 13 is given in Figures 16-18. In these figures the theo- 
retical curve is for d/h=0.5? the experimental results were 
obtained for 0. 492 < d/h < 0 . 511. The figures show_the comparison 
for three values of a, a = 0, a =0.05, a =0.10 in a-cjy(l+a). 

From Figures 14-18 it may be seen that the concept of plas- 
tic instability appears to be rather effective in correlating 
the results of the rupture tests. The figures also indicate tnat 
in the type of analysis which is based on an elastic-perfect y- 
plastic strip model, some increase in the yield strength (to a 
value called flow stress) is necessary. It should be pointed ^ 
out that even though the critical crack opening stretch criterion 
qave rather good results for aluminum and titanium cylinders 
(see Figure 12 and the work of Erdogan and Ratwani ) , it was not 
satisfactory for steel cylinders* ** . 


5 . Depressurization and Crack Arrest 


In pressurized pipes, containers, and other shell struc 
tures, generally the first phase of fracture. is the ru P t “^® 
the net ligament around a defect zone shown in Figure 2 through 
any one or a combination of fracture mechanisms.. This formation 
of a through crack will be followed by depressurization. 
ensuing phenomenon leads to a coupled fracture dynamics - fluid 
mechanics problem which is rather difficult to formulate. It is 
worthwhile, however, to mention two idealized cases for which 1 
may be possible to construct approximate models. 

The first refers to an infinitely long cylindrical shell 
containing pressurized gas. If the velocity of sound in Jf* s 
is smaller than the terminal velocity of the crack in the shell 
(which is usually the case) , and if after the formation of the 
through crack the initial acceleration phase of the crack. propa 
qation is of sufficiently short duration, then it is possible 
for the crack to grow at a constant velocity driven^by a constant 
pressure (which is the initial pressure) . In this steady-state 
case the elastic energy released by the shell and the expanding 
gas will be spent largely as kinetic energy (for the unfolding 
shell walls and for the gas) , plastic work, and fracture energy, 
the energy balance equation determining the crack velocity. 
Because of the geometry of the medium, the stress waves generated 
alonq the crack path in the shell will interfere with the mam 
membrane loads after traveling around the shell circumference. 
Hence, even in this- "steady-state" crack propagation a certain 
amount of waviness in the crack path would be expected. If the 
driving force is sufficiently high, it may not be possible to 
maintain a steady longitudinal crack path,^and the crack will 
keep branching to absorb the excess energy . On the other hand, 

* The explanation of the success or failure of the two criteria 
in their application to the two groups of materials under con- 
sideration may perhaps lie in the differences in the nature o 
the plastic deformations and in the shape of the yield zones 
around the crack tip. 

**The natural waviness of the crack path and the material imper- 
fections will also increase the tendency of the crack to branch 
or deviate from a meridional path. 
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if the driving force is not high enough, the crack velocity 
corresponding to the energy balance may^be less than ^ the soun 
velocity in the gas. In this case the steady state crack 
propagation cannot take place and as a result of depressuriza- 
tion the crack. would eventually be arrested. 

The second problem refers to the fracture of a finite 
cylindrical container. In this case one may consider the problem 
as being. quasi-static. The justification of this assumption is 
based on the fact that the part of the phenomenon which is of 
primary interest here is the initial, stage for which the cracK 
velocity and as a result the inertia effects are relative * , 
small, and the distortion in the overall geometry of ^ the < shell 
is negligible (meaning that the existing shell solution is still 
valid) . If this phase of the study leads to the conclusion th 
the fracture Will not be arrested and will continue with mcrea 
ing velocity and overall distortion, the inertia effects may 
indeed not be negligible. However, m this case the failure is 
very likely to be catastrophic and the solution, even if it can 
be obtained, is of little practical interest. 

The propagation and arrest of a crack in a finite volume 
container will now be considered under the following additions 
assumptions: 

(a) The medium in the container behaves either as a per- 
fect gas or as a nonviscous liquid. 

(b) The pressure drop in the container is assumed to be 
sufficiently "slow" so that the fluid, mechanics problem may be 
treated as "one-dimensional" with the basic variables p (pres- 
sure) , p (density) , and v (discharge velocity) . 

(c) The ’’orifice 1 ' may be approximated by an ellipse with 
the axes 2a (the crack length) and 6 0 (the COD at x-0). 

(d) The last assumption concerns the dynamic rupture cri- 
terion or some statement about the motion of the crack which 
provides one of the two basic equations for the solution of the 
problem. Very little is known quantitatively about the velocity 
dependence of the fracture resistance parameters such as K IC or 
S- aside from the fact that they are sensitive to the crack 
velocity as well as to the temperature. Therefore in the model 
that follows it will simply be assumed that after the rupture 
the net ligament the crack will grow either with a constant vel- 
ocity or with a constant acceleration. However, the analysis is 
such that one may postulate any function for the crack lengt , 

2a (t) . 

' * Within the confines of these assumptions, the basic vari- 

ables of the problem are p(t) (pressure), p (t) (density) ,v(t) 
(discharge velocity), and a(t) (the half-crack length). The 
relevant equations to determine these functions are: 

(i) pP~^ = constant (adiabatically expanding gas) , or 


PP 


= constant ( isothermally expanding liquid) , 
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(ii) The energy equation which may be expressed as (see 
Landau and Lifshitz 10 ) 

v = (_2l_ E) 1/2 (gas), v= [2(p-p„)/p„] 1/2 (liquid). 
Y-l p e e 

(10) 

(iii) The continuity equation 

xdp = - ftA(t) p e v(t)dt (11) 

where t is the volume of the container, ft is a discharge shape 
factor, A(t) is the area of the ’'orifice", and p Q is the density 
at the "exit"* From the shell solution A may be approximated by 


. , . % 2 , 2R . 1. 47a" * 

A{t) = ira p ( rrr + ^ ) - 


hE 


h 2 E 


(iv) Motion of the crack 

a(t) » a + v t (constant velocity), 
o c 


( 12 ) 


(13) 


a(t) 


a + i bt 2 (constant acceleration) , 
o 2 


or any given function. 


(14) 


Solving these equations for a constant velocity crack in a 
gas-filled container we find 


p(a) = { P :“ + 


7taC 

hEv 


2R, 3 3. . 1.47, 5 15. ,,-!/« 


[^(a J - ap + 


5h 


(15) 


„ Oy / 2 y ~ _l/2y„l/Y 

t v iy< 


■l)p (0T P o p e 


a - 3 (Y~l) 

a " 2y 


where p Q is the initial pressure, p Q = p(a Q ). 

Similarly for a constant acceleration crack we obtain 


p (a) = [p~“ + aKf(a)3~ 1/a , 


(16) 


K = 


1. 47 ttC 
h 2 E/2b 


b, = 


2Rh 


1 1.47 ' 


. _ 4a r- b-, + 6a « c 

f(a , = ^( a - a ^ 4 - 5 + 3 ^( a - a o ) 3 - 5 + ^-^( a - a ,) 2 * 5 


2a o (b l +2a o ) 

1.5 


(a ” a 0 > + 2 a 2 (b 1 + a 2 ) (a- a o ) 0,5 


Similarly for a container filled with a nonviscous liquid 
under isothermal conditions and for a constant velocity crack 
we obtain 
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p(a> - P = { (P n - Pj 


- 1/2 


/X [f(a 3 -a 3 ) 

EhTV c V p e 3 ° 


(a 5 - a^) 3 } 


An expression similar to (16) may be obtained for a constant 
acceleration crack. 

Figure 19 shows qualitatively the variation of pressure for 
a constant velocity crack (a = constant) , for a constant accelera- 
tion crack (a = constant), and that corresponding to a constant 
crack opening stretch 6 = 6 C as a functionof the crack length 
(see Figure 9). 6 =<5 C represents the resistance of./the shell 

whereas the depressurization curves a - constant or a const 
represent the external loads. When the depressurization curve 
falls below the resistance curve the crack will be arrested. The 
main feature of the depressurization curves is that the possibil 
ity of crack arrest increases with increasing <5 C and decreasing 
crack velocity v c or crack acceleration b. 

Figure 20 shows a sample quantitative example for a constant 
velocity crack. In the figure the solid line corresponds to the 
static resistance curve (6^) =2. P Y = Oyh/R is the Pressure 
corresponding to fully yielded shell. The net ligament may break 
through at a pressure lower than that which would initiate static 
fracture at the tip of a through crack of length 2a. The drop in 
the resistance curve (dashed line) shown in the figure is assumed 
to account for the dynamic effects. Hence, when the net ligament 
ruptures, the input pressure p Q will generally be somewhat . higher 
than the resistance (or equilibrium) value and the crack will 
start or continue to grow. The information used_ in the problem 
is: ay = 70,000 psi, p Q = 434 psi, a 0 = 1 in., *c “ 2d l' ? h " J? * 1 

r = 10 in., 1 = 20 in. (length of the container), E - 10 psi, 
v - 1 4 p- = 14.7 psi, v c = 5 in. /sec., ft-0.3 (for air), 0.4 
(for* water) . a = a r shown in the figure corresponds to the arrest 
value of the crack length in the liquid- filled container. The 
figure indicates that even for a very slow crack velocity assumed 
in the solution, for an expanding gas the crack arrest (before 
the fracture reaches catastrophic proportions) does not seem to 
be possible. However, for expanding liquids the crack arrest 
during the initial stage of crack propagation appears to be quite 

possible. 

In the foregoing example the critical crack opening stretch 
was used as the fracture criterion. It should however, be noted 
that the problem could be solved by using the concept of plastic 
instability or any other fracture criterion and by assuming a 
more general function a(t> for the length of the propagating 

crack. 
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Fig. 1. Part-through crack geometries. 



Fig. 2. Internal flaw in the shell with fully-yielded net 
ligaments . 






Fiq. 3. Crack opening displacement in the neutral surface of 
the shell and at the mid-section of the crack (x-0 r z-0). 


X=0.5 


6 


Fig. 4. Crack opening stretch in the neutral surface and at 
the crack tip (x=a r z=0). 






d/h 

Pig. 5. N 0 vs. d for constant total crack opening stretch, 
at the leading edge of the crack in the mid-section {x=0, 
z = h/2 - d) in a cylindrical shell with an external surface 
crack, A = 1. 


Fig. 6. N 0 vs. d for constant total crack opening stretch, 
at the leading edge of the crack in the mid-section {x=0, 
z = h/2 - d) in a cylindrical shell with an external surface 
crack, A = 2 . 
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Fig. 7. N 0 vs. d for constant total crack opening stretch, 6t 
at the leading edge of the crack in the mid-section (x— 0, 
z = h/2 - d) in a cylindrical shell with an external surface 
crack, X = 3 . 


Fig. 8. N 0 vs. d for constant total crack opening stretch, <5 
at the leading edge of the crack in the mid-section {x= 0, 
z=h/2-d) in a cylindrical shell with an external surface 
crack, A = 4. 
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Fig. 9. 
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pig. ii. Comparison of the calculated COD values with the 
experimental measurements of Cowan and Langford (1969) . 
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Fig. 13. Load carrying capacity of pressurized cylindrical 
shells with a part-through or through crack based on the 
plastic instability criterion. 
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Fig. 14. Comparison of the results of the fracture tests in 
steel pipes with a through crack (Kiefner, et al. 1972} with 
that given by the plastic instability criterion. 
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Fig. 15. Same 
Kiefner, et al. 

as Fig. 14 
(1972) . 

for 

all 

steel 

pipes tested by 


Thaor. {$=0.5) 
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Fig. 16. The results of the fracture testa in steel pipes with 
a part-through crack, 0.492 < d/h < 0.511, o^l-loy {Kiefner, et 
al. 1972). Solid line is the result predicted by fracture 
instability for d/h =0.5. 
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Fig. 19. Qualitative nature of the fracture resistance (« - ' , 

and constant velocity and constant acceleration depressurizatio 
curves. 2a Q : initial length of the surface flaw; p Q . initial 
pressure. 
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Fig. 20. Pressure vs. crack length corresponding to a constant 
crack opening stretch 6 =6 C (= 2d].) at the crack tip (the resist- 
ance curves) , and liquid and gas depressurization curves in a 
finite-volume cylindrical container with a propagating through 
crack . 



